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his paper provides a review of the variants of dark matter (CDM, HDM) which are 
thought to be fundamental components of the universe and their role in origin and 
evolution of structures. 



^ ' Introduction 

oo : 

The origin and evolution of large scale structure is today the outstanding problem in 
' cosmology. This is the most fundamental question we can ask about the universe whose 
solution should help us to better understand problems as the epoch of galaxy formation, 
the clustering in the galaxy distribution, the amplitude and form of anisotropics in the 
^ microwave background radiation. Several has been the approaches and models trying to 

\ attack and solve this problem: no one has given a final answer. 
C3 The leading idea of all structure formation theories is that structures was born from small 
^ • perturbations in the otherwise uniform distribution of matter in the early Universe, which 
' is supposed to be, in great part, dark (matter not detectable through light emission). 
With the term Dark Matter cosmologists indicate an hypothetic material component 
of the universe which does not emit directly electromagnetic radiation (unless it decays 
in particles having this property (Sciama 1989, but also see Bowyer et al. 1999)). 
Dark matter, cannot be revealed directly, but nevertheless it is necessary to postulate 
its existence in order to explain the discrepancies between the observed dynamical 
properties of galaxies and clusters of galaxies and the theoretical predictions based upon 
models of these objects assuming that the only matter present is the visible one. If in the 
space were present a diffused material component having gravitational mass, but unable 
to emit electromagnetic radiation in significative quantity, this discrepancy could be 
eliminated (Turner 1991). The study of Dark Matter has as its finality the explanation 
of formation of galaxies and in general of cosmic structures. For this reason, in the last 
decades, the origin of cosmic structures has been "framed" in models in which Dark 
Matter constitutes the skeleton of cosmic structures and supply the most part of the 
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mass of which the same is made. 

There are essentially two ways in which matter in the universe can be revealed: by 
means of radiation, by itself emitted, or by means of its gravitational interaction with 
baryonic matter which gives rise to cosmic structures. Electromagnetic radiation permits 
to reveal only baryonic matter. In the second case, we can only tell that we are in 
presence of matter that interacts by means of gravitation with the luminous mass in 
the universe. The original hypotheses on Dark Matter go back to measures performed 
by Oort (1932) of the surface density of matter in the galactic disk, which was obtained 
through the study of the stars motion in direction orthogonal to the galactic plane. 
The result obtained by Oort, which was after him named "Oort Limit", gave a value of 
p = 0.15M o pc" 3 for the mass density, and a mass, in the region studied, superior to that 
present in stars. Nowadays, we know that the quoted discrepancy is due to the presence 
of HI in the solar neighborhood. Other studies (Zwicky 1933; Smith 1936) showed the 
existence of a noteworthy discrepancy between the virial mass of clusters (e.g. Coma 
Cluster) and the total mass contained in galaxies of the same clusters. These and other 
researches from the thirties to now, have confirmed that a great part of the mass in the 
universe does not emit radiation that can be directly observed. 



1.1 Determination of SI and Dark Matter 



The simplest cosmological model that describes, in a sufficient coherent manner, the 
evolution of the universe, from 10 -2 s after the initial singularity to now, is the so called 
Standard Cosmological Model (or Hot Big Bang model). It is based upon the Friedmann- 
Robertson- Walker (FRW) metric, which is given by: 
dr 2 



ds 2 = c z dt z - a(tf 



+ r z (d0 z + sinO z d<j) Z 



(1) 
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where c is the light velocity, a(t) a function of time, or a scale factor called "expansion 
parameter" , t is the time coordinate, r, 9 and <fr the comoving space coordinates. The 
evolution of the universe is described by the parameter a(t) and it is fundamentally 
connected to the value p of the average density. 

The equations that describe the dynamics of the universe are the Fricdmann's equations 
(Friedman, 1924) that we are going to introduce in a while. These equations can be 
obtained starting from the equations of the gravitational field of Einstein (Einstein, A., 
1915): 
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R-ik — -QikR = —Tik (2) 

A C 



where now, i?^ is a symmetric tensor, also known as Ricci tensor, which describes the 
geometric properties of space-time, gik is the metric tensor, R is the scalar curvature, Tik 
is the energy-momentum tensor. 



These equations connect the properties of space-time to the mass-energy. In other terms 
they describe how space-time is modeled by mass. Combining Einstein equations to the 
FRW metric leads to the dynamic equations for the expansion parameter, a(t). These 
last are the Friedmann equations: 

d(pa 3 ) = -pd(a 3 ) (3) 
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2^ + % + 4 = -SttGp (5) 

where p is the pressure of the fluid of which the universe is constituted, k is the 
curvature parameter and a(t) is the scale factor connecting proper distances r to the 
comoving ones x through the relation r = a(<)x. One of the components of the today 
universe are galaxies. If we assume that galaxies motion satisfy Weyl (Weyl, H., 1923) 
postulate, the velocity vector of a galaxy is given by u l — (1,0,0,0), and then the 
system behaves as a system made of dust for which we have p — 0. Only two of the 
three Friedmann equations are independent, because the first connects density, p to the 
expansion parameter a(t). The character of the solutions of these equations depends 
on the value of the curvature parameter, k, which is also determined by the initial 
conditions by means of Eq. 3. The solution to the equations now written shows that if 
p is larger than p c = = 1.88 * 10~ 29 g / cm 3 (critical density, which can be obtained 
from Friedmann equations putting t = to, k = 0, and H = lOOfcm/ sMpc), space-time 
has a closed structure (k — 1) and equations shows that the system go through a 
singularity in a finite time. This means that the universe has an expansion phase until it 
reaches a maximum expansion after which it recollapse. If p < p c , the expansion never 
stops and the universe is open k = — 1 (the universe has a structure similar to that of an 
hyperboloid, in the two-dimensional case). If finally, p = p c the expansion is decelerated 
and has infinite duration in time, k = 0, and the universe is flat (as a plane in the 
two-dimensional case). The concept discussed can be expressed using the parameter 
Q = In this case, the condition SI = 1 corresponds to k = 0, fi < 1 corresponds to 

Pc 

k = — 1, and SI > 1 corresponds to k = 1. F] 

The value of SI can be calculated in several ways. The most common methods are the 
dynamical methods, in which the effects of gravity are used, and kinematics methods 
sensible to the evolution of the scale factor and to the space-time geometry. The results 
obtained for SI with these different methods are summarized in the following. 



Dynamical methods: 



See next paragraphs for some items on cosmological models with non-zero cosmological 
constant. 



(a) Rotation curves: The contribution of spiral galaxies to the density in the universe 
is calculated by using their rotation curves and the third Kepler law. Using the last it is 
possible to obtain the mass of a spiral galaxy from the equation: 

M(r) = v 2 r/G (6) 

where v is the velocity of a test particle at a distance r from the center and M(r) is 
the mass internal to the circular orbit of the particle. In order to determine the mass 
M is necessary to have knowledge of the term v 2 in Eq. (§) and this can be done 
from the study of the rotation curves through the 21 cm line of HI. Rotation curves of 
galaxies are characterized by a peak reached at distances of some Kpcs and a behavior 
typically flat for the regions at distance larger than that of the peak. A peculiarity is 
that the expected Keplerian fall is not observed. This result is consistent with extended 
haloes containing masses till 10 times the galactic mass observed in the optical (van 
Albada, T. S., Sancisi, 1986). The previous result is obtained assuming that the halo 
mass obtained with this method is distributed in a spherical region so that we can use 
Eq. (^) and that we neglect the tidal interaction with the neighboring galaxies which 
tend to produce an expansion of the halo. After M and the luminosity of a series of 
elliptical galaxies is determined, the contribution to the density of the universe is given 
by p =< > I where I is the luminosity per unit volume due to galaxies and can be 
obtained from the galactic luminosity function (f>(L)dL, which describes the number of 
galaxies per Mpc 3 and luminosity range L, L + dL. The value that is usually assumed for 
I is I = 2Ahl0 8 Li, o Mpc~ 3 . The arguments used lead to a value of fl g for the luminous 
parts of spiral galaxies of fl g < 0.01, while for haloes flh > 0.03 — 0.1. The result shows 
that the halo mass is noteworthy larger than the galactic mass observable in the optical 
(Peebles,P.J.E.,1971). 

(b) Virial theorem: 

In the case of non spiral galaxies and clusters, the mass can be obtained using the 
virial theorem 2T + V = 0, with 

T = Im < v 2 > (7) 

where < v 2 > is the velocity dispersion along the line of sight. After getting the value 
of M of the cluster by means of the virial theorem one determines L by means of 
observations. Given M and L, the value of O for clusters is obtained similarly to the 
case of spiral galaxies. Usual values obtained for f2 are = 0.1 — 0.3 (Peebles, P.J.E., 
1971). A problem of the quoted method is that in general the results obtained are the 
right one only for virialized, spherically symmetric clusters. In general, clusters are not 
virialized objects: even Coma clusters seems to have a central core constituted by more 



than one blob of mass (Henry, J.P.,Briel, U.G.,1992). 
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(c) Peculiar velocities: 



The velocity of a galaxy can be written as: 

V g = Hr + <5v (8) 

where H is Hubble constant. The previous equation shows that the motion of a galaxy 
is constituted by two components: the velocity of the galaxy due to the Hubble flow 
and a peculiar velocity <$v, which describes the motion of the galaxy with respect to the 
background. In the linear regime, as we see in a while, we find that on average on a scale 
of length A, it is: 

— W £1°' 6 — ^-y— (9) 

c H^ 1 p 

(Kolb c Turner 1990). Then given the overdensity on scale A and 5v, it is possible to 
obtain Q. The overdensity ^ can be obtained from the overdensity of galaxies using 
the relation — = -^-b^ 1 with 1 < b < 3. Using IRAS catalog in order to obtain the 
overdensity in galaxies one finds Q = 1 ( M.S. Turner ). The values of Q obtained using 
the method of peculiar velocity assume that the peculiar velocity fields describe in an 
accurate way, the inhomogeneity in the distribution of underlying mass. We should note 
that the peculiar velocity method has some difficulties. In general, in order to obtain 
these last it is necessary to determine the redshift and the distance of galaxies and by 
using these data it is possible to obtain the peculiar velocity: 



u pec 



zc~H d (10) 



It is evident that there are problems in measuring the distance d, problems connected 
to difficulties in finding trustable indicators of distance. Moreover the peculiar velocity 
can be determined only along the line of sight. 



(d) Kinematic methods: 



These methods are based upon the use of relations between physical quantities depen- 
dent on cosmological parameters. An example of those relations is the relation luminosity 
distance-rcdshift : 

H d L = z + ^(l-q )z 2 (11) 

where H is Hubble constant nowadays, z is the redshift, d^ = ^r- the luminos- 
ity distance, L the absolute luminosity, and F the flux. By means of the relations 
luminosity-redshift, angle- redshift, number of objects-redshift, it is possible to determine 
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the parameter of deceleration qo = — rfi aQ (do an d Hq = = lQOhkm/Mpcs arc the 
scale factor and the Hubble constant, nowadays). At the same time qo is connected to 
Q by means of go = §j m a matter universe. One of the first test used, the luminosity 
distance-redshift has several problems due to effects of the evolution of sources. 
Uncertainties in the knowledge of the effects of galactic evolution on the intrinsic 
luminosity of the same has not, in the past, permitted to find definitive values of qo. 
For this reason, Loh & Spillar (1986) introduced another kinematic test: number of 
galaxies-redshift. This test is based on the count of the number of galaxies in a comoving 
element of galaxies, defined by the surface dQ and the redshift dz. This number depends 
on qo. Nevertheless the effects of evolution of sources influences on the results of the 
test, it is more sensible to the evolution of number of sources than to the evolution of 
luminosity, on which there is not an accepted theory. Results gives high values of 
(fi = Q.9to 7 5 ) (Kolb e Turner 1990). 

(e) Primordial nucleosynthesis: 
The theory of primordial nucleosynthesis, proposed in 1946 by Gamov, assumes that the 
light elements till Li-j are generated after big bang and that heavier elements originate 
from nuclear reactions inside stars. The values obtained for the abundances depends 
on some parameters like: ?y, the value of the ratio baryons-photons, nowadays; N v , the 
number of neutrinos species; Tcmbr-, the temperature of Cosmic Microwave Background 
Radiation. The theory of primordial nucleosynthesis permits to give limits to fib (b 
stands for baryons). With a ratio baryons-photons 3 * 10 -10 < T) < 5 * 1CP 10 , and a value 
of N v < 4 for the neutrinos species, T C mbr = 2.736 ±0.01AT, is found 0.011 < tt b < 0.12 
(Kolb e Turner 1990 ). 

In the following, we summarize the results of some more recent results. 

Dekel et al 1996 used several methods to obtain the value of SI. According to their 
classification, we divide the methods into the following four classes: 

• Global measures. Based on properties of space-time that constrain combinations of 
£l m and the other cosmological parameters (A, Ho, to). 

• Virialized Systems. Methods based on nonlinear dynamics within galaxies and 
clusters on comoving scales 1 — 10h~ 1 Mpc. 

• Large-scale structure. Measurements based on mildly-nonlinear gravitational dy- 
namics of fluctuations on scales 10 — 100h~ 1 Mpc of superclusters and voids, in 
particular cosmic flows. 



2 Among kinematics methods we should mention SNela which played a key role in the last 
few years (Valageas 1999). 
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• Growth rate of fluctuations. Comparisons of present day structure with fluctuations 
at the last scattering of the cosmic microwave background (CMB) or with high 
redshift objects of the young universe. 
The methods and current estimates are summarized in Table 1. The estimates based 
on virialized objects typically yield low values of fl m ~ 0.2 — 0.3. The global measures, 
large-scale structure and cosmic flows typically indicate higher values of Q m ~ 0.4 — 1. 
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Table 1: Estimates of fi 
Table 1: Estimates of VL m 



Global Measures 



Inflation, Occam 
Lumimosity distance SNIa 

Lens Counts 
CMB Peak 



n m + n A = i (n m = i, n A = o) 

-0.3 < a m - fl A < 2.5 (90%) 
Flat Q m > 0.49 (95%) 
Flat Q m > 0.34 (95%) 

n m + n A < 1.5 (95%) 

n m + fl A > 0.3 (likely - 0.7) 
Q m - 0.7n A < 1.3 (likely < 0) 



Virialized Objects 



(M/L)C 
Baryon fraction 
Cosmic Virial Th. 
Local Group 



n m ~ 0.25 (0.1 - 1.0) 
On^s 2 ~ 0.3 - 0.5 (low-high fi 6 ) 
Point mass £l m ^0.2 (halos— >1) 
Point mass fi m ~0.15 (halos — >0.7) 



Large-Scale; Flows Peculiar velocities 

Rcdshift Distortions 
Velocity vs Density 

Cluster Abundance 



il m > 0.3 (95%) 

Q^Vf = 0.8 ± 0.2 (/3 b j ~ 1.05 c ) 
/?/ - 0.5-1.2 

0.5 — 1.2 (scale dependent) 
/3 ~ 0.4-0.95 

tt°Ja B ~ 0.5 - 0.6 (/?/ ~ 0.7- 0.8 C ) 



Fluctuation Growth Cluster Morphology 
Galaxy Formation 

P fe (p) VS C; 

P fe («) vs C, 



« m > 0.2 (?) 

(?) 

CDM n = 1 b = 1: O m /i 6 5 - 0.3 
CDM flat: n m h 65 n 2 ~ 0.7 ± 0.1 



° See Dekel et al. 1996 for a discussion of the table. <r 8 is the rms mass density fluctuation 
in a top- hat sphere of radius 8h~ 1 Mpc. 

b [3 = fl°- 6 /b, bj for IRAS galaxies, b Q for optical galaxies. 

c b /b!~ 1.3, b ~l/a 8 . 

Neta Bahcall et al. 1997, showed that the evolution of the number density of rich 
clusters of galaxies breaks the degeneracy between Q (the mass density ratio of the 
universe) and a 8 (the normalization of the power spectrum), as ^°' 5 — 0.5, that follows 
from the observed present-day abundance of rich clusters. The evolution of high-mass 
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(Coma-like) clusters is strong in 51 = 1, low-ag models (such as the standard biased 
CDM model with erg — 0.5), where the number density of clusters decreases by a factor 
of ~ 10 3 from z = to z ~ 0.5; the same clusters show only mild evolution in low-fi, 
high-ag models, where the decrease is a factor of ~ 10. This diagnostic provides a most 
powerful constraint on O. Using observations of clusters to z ~ 0.5 — 1, the authors 
found only mild evolution in the observed cluster abundance, and Q = 0.3 ±0.1 and 
(T 8 = 0.85 ± 0.15 (for A = models; for Cl + A = 1 models, Cl = 0.34 ± 0.13). 

Ferreira et al. 1999, proposed an alternative method to estimate V12 directly from 
peculiar velocity samples, which contain redshift-indcpcndcnt distances as well as galaxy 
rcdshifts. In contrast to other dynamical measures which determine /3 = f2°' 6 er 8 , this 
method can provide an estimate of O 06 cr| for a range of a% where Q is the cosmological 
density parameter, while as is the standard normalization for the power spectrum of 
density fluctuations. 

Melchiorri et al. 1999, used the angular power spectrum of the Cosmic Microwave 
Background, measured during the North American test flight of the BOOMERANG 
experiment, to constrain the geometry of the universe. Within the class of Cold Dark 
Matter models, they find that the overall fractional energy density of the universe, fi, is 
constrained to be 0.85 < fi < 1.25 at the 68% confidence level. 

Branchini et al. 1999, compared the density and velocity fields as extracted from the 
Abell/ACO clusters to the corresponding fields recovered by the POTENT method from 
the Mark III peculiar velocities of galaxies. Quantitative comparisons within a volume 
containing <~ 12 independent samples yield (3 C = Sl°- 6 /fe c = 0.22 ± 0.08, where b c is the 
cluster biasing parameter at 15h~ 1 Mpc. If b c ~ 4.5, as indicated by the cluster correlation 
function, their result is consistent with Q <~ 1. 

(f) Inflation: 

It is widely supposed that the very early universe experienced an era of inflation (see 
Guth 1981, Linde 1990, Kolb & Turner 1990). By 'inflation' one means that the scale 
factor has positive acceleration, a > 0, corresponding to repulsive gravity and 3p < —p. 
During inflation aH = a is increasing, so that comoving scales are leaving the horizon 
(Hubble distance) rather than entering it, and it is supposed that at the beginning of 
inflation the observable universe was well within the horizon. 

The inflationary hypothesis is attractive because it holds out the possibility of cal- 
culating cosmological quantities, given the Lagrangian describing the fundamental inter- 
actions. The Standard Model, describing the interactions up to energies of order ITeV, 
is not viable in this context because it does not permit inflation, but this should not be 
regarded as a serious setback because it is universally agreed that the Standard Model 
will require modification at higher energy scales, for reasons that have nothing to do with 
cosmology. The nature of the required extension is not yet known, though it is conceiv- 
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able that it could become known in the foreseeable future (Weinberg 1993). But even 
without a specific model of the interactions (ie., a specific Lagrangian), the inflationary 
hypothesis can still offer guidance about what to expect in cosmology. More dramatically, 
one can turn around the theory-to-observation sequence, to rule out otherwise reasonable 
models. The importance of inflation is connected to: 

a) the origin of density perturbations, which could originate during inflation as quantum 
fluctuations, which become classical as they leave the horizon and remain so on re-entry. 
The original quantum fluctuations are of exactly the same type as those of the electro- 
magnetic field, which give rise to the experimentally observed Casimir effect. 

b) One of the most dramatic and simple effects is that there is no fine-tuning of the 
initial value of the density parameter fl = Sirp/Sm^H 2 . From the Friedmann equation, 
Q, is given by 

ft-l = (J-:) 2 (12) 

Its present value Qq is certainly within an order of magnitude of 1 , and in the absence of 
an inflationary era fl becomes ever smaller as one goes back in time, implying an initial 
fine tuning. In contrast, if there is an inflationary era beginning when the observable 
universe is within the horizon, Eq. (|l^) implies that Slo will be of order 1, provided only 
that the same is true of at the beginning of inflation. A value of f2o extremely close to 
1 is the most natural, though it is not mandatory.^] 

c) Another effect of inflation is that it can eliminate particles and topological defects 
which would otherwise be present. Anything produced before inflation is diluted away, 
and after inflation there is a maximum temperature (the 'reheat' temperature) which is 
not high enough to produce all the particles and defects that might otherwise be present. 
As we shall remark later, this mechanism can remove desirable, as well as undesirable, 
objects. 

d) The most dramatic effect of inflation is that it may offer a way of understanding the 
homogeneity and isotropy of the universe, or at any rate of significant regions of it. We 
have nothing to say about this complex issue in its full generality, but a more modest 
version of it is our central concern. In this version, one begins the discussion at some early 
stage of inflation, when the universe is supposed already to be approximately homoge- 
neous and isotropic. One then argues that in that case, scales far inside the horizon must 
be absolutely homogeneous and isotropic, except for the effect of vacuum fluctuations in 
the fields. Finally, one shows that after they leave the horizon, such fluctuations can be- 
come the classical perturbations that one deals with in cosmological perturbation theory. 
This possibility was first pointed out for gravitational waves by Starobinsky (1979) and 



3 An argument has been given for fio very close to 1 on the basis of effects on the cmb 
anisotropy from regions far outside the observable universe (Turner 1991), but it is not valid as 
it stands because it ignores spatial curvature. 
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for density perturbations by several people (Guth & Pi 1982; Hawking 1982; Starobinsky 
1982). As we shall go to some trouble to demonstrate, the vacuum fluctuations can be 
evaluated unambiguously once an inflationary model is specified, 
(g) Scalar field inflation: 

Two mechanisms for inflation have been proposed. The simplest one (Guth 1981) 
invokes a scalar field, termed the inflaton field. An alternative (Starobinsky 1980) is to 
invoke a modification of Einstein gravity, and combinations of the two mechanisms have 
also been proposed. During inflation however, the proposed modifications of gravity can 
be abolished by redefining the spacetime metric tensor, so that one recovers the scalar 
field case. We focus on it for the moment, but modified gravity models will be included 
later in our survey of specific models. 

In comoving coordinates a homogeneous scalar field <f> with minimal coupling to grav- 
ity has the equation of motion 

4> + 3H4> + V'(4>) = (13) 
Its energy density and pressure are 

P = v+ l -4> 2 (14) 
P=-V+\tf (15) 

If such a field dominates p and p, the inflationary condition 3p < p is achieved provided 
that the field rolls sufficiently slowly, 

<p 2 < V (16) 

Practically all of the usually considered models of inflation satisfy three conditions. 
First, the motion of the field is overdamped, so that the 'force' V balances the 'friction 
term' 3H<fr, 

**-±r <i7) 



Second, 



1 pi 



V , 

77 «1 (18) 



16tt V V , 

which means that the inflationary requirement 4> 2 < V is well satisfied and 

H 2 - \^rV (19) 

These two conditions imply that H is slowly varying, and that the scale factor increases 
more or less exponentially, 

a oc e Ht (20) 
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The third condition that is usually satisfied is 

M « 1 (21) 
where 

"=-£ir (22) 



It can be 'derived' from the other two by differentiating the approximation Eq. (|17|) for 
4> and noting that consistency with the exact expression Eq. ([l3]) requires if) V is 
satisfied. However there is no logical necessity for the derivative of an approximation 
to be itself a valid approximation, so this third condition is logically independent of 
the others. Conditions involving higher derivatives of V could be 'derived' by further 
differentiation, with the same caveat, but the two that we have given, involving only the 
first and second derivatives, are the ones needed to obtain the usual predictions about 
inflationary perturbations. The term 'slow-roll inflation' is generally taken to denote a 
model in which they are satisfied and we are adopting that nomenclature here. Practically 
all of the usually considered models of inflation satisfy the slow-roll conditions more or 
less well. 

It should be noted that the first slow-roll condition is on a quite different footing from 
the other two, being a statement about the solution of the field equation as opposed to a 
statement about the potential that defines this equation. What we are saying is that in 
the usually considered models one can show that the first condition is an attractor solu- 
tion, in a regime typically characterized by the other two conditions, and that moreover 
reasonable initial conditions on (f> will ensure that this solution is achieved well before 
the observable universe leaves the horizon. 

(h) Conclusions: 

We have seen the possible values of ft using different methods. We have to add that 
Cosmologists are "attracted" by a value of fig = 1. This value of is requested by 
inflationary theory. The previous data lead us to the following hypotheses: 

i) Qo < 0.12; in this case one can suppose that the universe is fundamentally made of 
baryonic matter (black holes; Jupiters; white dwarfs). 

ii) f2o > 0.12; in this case in order to have a flat universe, it is necessary a non-baryonic 
component, fib = 1 is excluded by several reasons (see Efstathiou 1990 Kolb e Turner 
1990). The remaining possibilities are: 

1) existence of a smooth component with fl = 0.8. 

The test of a smooth component can be done with kinematic methods. 

2) Existence of a cosmological term, absolutely smooth to whom correspond an energy 
density p vac = 

3) existence of non-baryonic matter: the universe is fundamentally done of particles 
(neutrinos, WIMPS (Weakly Interacting Massive Particles)). 
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4) A combination of 2) and 3). 

Before going one, I want to recall that some authors (Finzi 1963; Milgrom 1983; Sanders 
1984 ) have assumed that we have a scant knowledge of physical laws. Sanders assumes 
that the gravitational potential changes with distance and in particular the gravitational 
constant has a different value at large distances. Milgrom assumes that the Newton 
law of gravitation is not valid when the gradient of the potential is small. In this case, 
the problem of the dynamics of clusters of galaxies is solved without introducing Dark 
Matter. In any case, the quoted assumptions have no general theory that can justify them. 

1.2 Dark matter in particles. 

We know that if f2 = 1, dark matter cannot be constituted exclusive of baryonic mat- 
ter. The most widespread hypothesis is that dark matter is in form of particles. Several 
candidates exist: neutrinos, axions, neutralinos, photinos, gravitinos, etc. Interesting par- 
ticles are usually grouped into three families: 

HDM (Hot Dark Matter), CDM (Cold Dark Matter) and WDM (Warm Dark Matter). 
In order to understand this classification it is necessary to go back to the early phases 
of universe evolution. The history of the universe is characterized by long phases of 
local thermodynamic equilibrium (LTE) and by "deviation" by it: nucleosynthesys, bari- 
ogenesys, decoupling of species, etc. In the early universe, were present the particles that 
we know today and other particles predicted theoretically, but that have not been ob- 
served. Massive particles preserved the thermodynamic equilibrium concentration until 
the rate, T, of reactions and interactions that produced that concentration was larger 
than the expansion rate of the universe, H. When the condition T > H was no longer 
satisfied the reactions stopped and the abundance of the considered species remained 
constant at the value it had at time of freez-out, Tf, time at which Y = H. If we indicate 
with Y = — the number of particles per unit comoving volume and we remember that n 
is the number density of species and s the entropy density, we obtain a contribution of 
the species to the actual density of the universe as flh 2 = 0.28Y(Tf)(—) (Kolb e Turner 
1990). At time Tf, particles could be relativistic or non-relativistic. Relativistic particles 
are today indicated with the term hot cosmic relics, HDM, while non-relativistic particles 
are named cold cosmic relics, CDM. There is an intermediate case, that of warm relics, 
WDM. 

An example of HDM are massive neutrinos. Possible masses for these neutrinos are: 



25ev < m v < lOOev 



(23) 



( Sciama, D.W., 1984) and 




4.9 — l.SGev for Maiorana's neutrinos 



1.3 — 4.2Gev for Dirac's neutrinos 



(24) 
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(Lee e Weinberg 1977 ). 

There are confirmed experimental evidence of the existence of massive neutrinos. In 
1980, Lyubimov et al announced the detection of an electronic antincutrino with mass 
30 eV, by means of the shape of the electron energy spectrum in the [3 decay of tritium 
(N.A.Jelley 1986 ). Experiments (Super-Kamiokande, SNOW) have obtained some 
evidence of non-zero mass from neutrino oscillations. This yields a difference of square 
masses of order 10~ 3 eV, and a mass of 0.05 eV (in the simplest case) (see Zuber 1998). 
Among typical examples of CDM we have WIMPS and in particular axions and 
neutralinos (SUSY particle). This particle was postulated in order to solve the strong 
CP problem in nuclear physics. This problem arises from the fact that some interactions 
violate the parity, P, time inversion, T, and CP. If these are not eliminated, they give rise 
to a dipole momentum for the neutron which is in excess of ten order of magnitude with 
respect to experimental limits (Kolb e Turner 1990). The solution to the problem was 
proposed by Peccei-Quinn in 1977 in terms of a spontaneous symmetry breaking scheme. 
To this symmetry breaking should be associated (as shown by Weiberg & Wilczec 1978) 
a Nambu-Goldstone boson: the axion. The axion mass ranges between 10~ 12 ev-lGev. In 
cosmology there are two ranges of interest: 10~ 6 ev < m a < 10 _3 ev ; 2>ev < m a < 8ev. 
Axion production in the quoted range can originate due to a series of astrophysical 
processes (Kolb & Turner 1990) and several are the ways these particles can be de- 
tected. Nevertheless the effort of researchers expecially in USA, Japan and Italy, axions 
remain hypothetical particles. They are in any case the most important CDM candidates. 

In the following, I am going to speak about the basic ideas of structure formation. I 
shall write about density perturbations, their spectrum and evolution, about correlation 
functions and their time evolution, etc. 

1.3 Origin of structures. 

Observing our universe, we notice a clear evidence of inhomogeneity when we consider 
small scales (Mpcs). In clusters density reaches values of 10 3 times larger than the average 
density, and in galaxies it has values 10 5 larger than the average density (Kolb & Turner 
1990). If we consider scales larger than 10 2 Mpcs universe appears isoptric as it is observed 
in the radio-galaxies counts, in CMBR, in the X background (Peebles 1971). The isotropy 
at the decoupling time, td ec , at which matter and radiation decoupled, universe was very 
homogeneous, as showed by the simple relation: 
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(Kolb e Turner 1990) []. The difference between the actual universe and that at 
decoupling is evident. The transformation between a highly homogeneous universe, at 
early times, to an highly local non homogeneous one, can be explained supposing that 
at tdec were present small inhomogeneities which grow up because of the gravitational 
instability mechanism (Jeans, J. H., 1902). Events leading to structure formation can be 
enumerated as follows: 

(a) Origin of quantum fluctuations at Planck epoch. 

(b) Fluctuations enter the horizon and they grow linearly till recombination. 

(c) Perturbations grow up in a different way for HDM and CDM in the post- 
recombination phase, till they reach the non-linear phase. 

(d) Collapse and structure formation. 

Before tdec inhomogeneities in baryonic components could not grow because photons and 
baryons were strictly coupled. This problem was not present for the CDM component. 
Then CDM perturbations started to grow up before those in the baryonic component 
when universe was matter dominated. The epoch t eq w 4.4 * 10 10 (£loh 2 )~ 2 sec, at which 
matter and radiation density are almost equal, can be considered as the epoch at which 
structures started to form. The study of structure formation is fundamentally an initial 
value problem. Data necessary for starting this study are: 

1) Value of Slo- In CDM models the value chosen for this parameter is 1, in conformity 
with inflationary theory predictions. 

2) The values of fli for the different components in the universe. For example in the case 
of baryons, nucleosynthesis gives us the limit 0.014 < < 0.15 while SIwimps ~ 0.9. 

3) The perturbation spectrum and the nature of perturbations (adiabatic or isocurva- 
ture). The spectrum generally used is that of Harrison- Zeldovich: P(k) = Ak n with 
n = 1. The perturbation more used are adiabatic or curvature. This choice is dictated 
from the comparison between theory and observations of CMBR anisotropy. 

1.4 The spectrum of density perturbation. 

In order to study the distribution of matter density in the universe it is generally 
assumed that this distribution is given by the superposition of plane waves independently 
evolving, at least until they are in the linear regime (this means till the overdensity 
S = << 1). Let we divide universe in cells of volume V u and let we impose periodic 
conditions on the surfaces. If we indicate with p the average density in the volume and 
with p(r) the density in r, it is possible to define the density contrast as: 

,( r) = (26) 



4 In fact, COBE data gives < 1(T S 
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This quantity can be developed in Fourier series: 

5(r) — Skexp(i'kr) = <5kea;p(— zkr) (27) 

k k 

( Kolb e Turner 1990), where k x = 2?r " x (and similar conditions for the other components) 
and for the periodicity condition 5(x,y,L) — 6(x 7 y,0) (and similar conditions for the 
other components). Fourier coefficients <5k are complex quantities given by: 

S u = — I 6(r)exp(-i\sr)dr (28) 

For mass conservation in V u we have also <5k=o = while for reality of <5(r), (5£ = i5_k- If 
we consider n volumes, V u , we have the problem of determining the distribution of Fourier 
coefficients 8k and that of \S\. We know that the coefficients are complex quantities and 
then 5k = |#k| exp(i8k). If we suppose that phases are random, in the limit V u — > oo it 
is possible to show that we get |<5| 2 = ^ k |<5k| 2 . The Central limit theorem leads us to 
conclude that the distribution for S is Gaussian: 
-S 2 

P(S) oc exp{^) (29) 



(Efstathiou 1990). The quantity a that is present in Eq. (29) is the variance of the density 
field and is defined as: 

^ 2 =<<5 2 >=£<N 2 >=^£^ ( 3 °) 
k Vu k 

This quantity characterizes the amplitude of the inhomogeneity of the density field. If 
V u — > oo, we obtain the more usual relation: 

° 2 = (2^f / mSk = i / p{k)k2dk (31) 

The term P(k) =< \S\ 2 > is called "Spectrum of perturbations". It is function only of k 
because the ensemble average in an isotropic universe depends only on r. A choice often 
made for the primordial spectrum is P(k) — Ak n which in the case n — 1 gives the 
scale invariant spectrum of Harrison-Zeldovich. An important quantity connected with 
the spectrum is the two-points correlation function £(r, t). It can be defined as the joint 
probability of finding an overdensity 5 in two distinct points of space: 

£(r,t) =< S(r,t)6(r + x,t) > (32) 

(Peebles 1980), where averages are averages on an ensemble obtained from several re- 
alizations of universe. Correlation function can be expressed as the joint probability of 
finding a galaxy in a volume 8V% and another in a volume 8V2 separated by a distance 

8 2 P = n 2 v [l+Z{r 12 )]SV 1 8V 2 (33) 

where ny is the average number of galaxies per unit volume. The concept of correlation 
function, given in this terms, can be enlarged to the case of three or more points. 
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Correlation functions have a fundamental role in the study of clustering of matter. If 
we want to use this function for a complete description of clustering, one needs to know 
the correlation functions of order larger than two (Fry 1982). By means of correlation 
functions it is possible to study the evolution of clustering. The correlation functions are, 
in fact, connected one another by means of an infinite system of equations obtained from 
moments of Boltzmann equation which constitutes the BBGKY (Bogolyubov-Green- 
Kirkwood-Yvon) hierarchy (Davis e Peebles 1978). This hierarchy can be transformed 
into a closed system of equation using closure conditions. Solving the system one gets 
information on correlation functions. 

In order to show the relation between perturbation spectrum and two-points correlation 
function, we introduce in Eq. (|32j), Eq. (|27|), recalling that <5£ = S(— k) and taking the 
limit V u — > oo, the average in the Eq. ( |32] ) can be expressed in terms of the integral: 

£(r) - 7^53 J \6(k)\ 2 exp(-ikr)d 3 k (34) 

This result shows that the two-point correlation function is the Fourier transform of the 
spectrum. In an isotropic universe, it is |r| = r and then |k| = k and the spectrum can 
be obtained from an integral on |k| = k. Then correlation function may be written as: 

«0 = i / k'Pikf-^dk (35) 

During the evolution of the universe and after perturbations enter the horizon, the spec- 
trum is subject to modulations because of physical processes characteristic of the model 
itself (Silk damping (Silk 1968) for acollisional components, free streaming for collisional 
particles, etc.). These effects are taken into account by means of the transfer function 
T(k]t) which connects the primordial spectrum P(k;t p ) at time t p to the final time tf\ 

P(k;t f ) = [^|] T 2 (k;t f )P(k;t p ) (36) 

where b(t) is the law of grow of perturbations, in the linear regime. In the case of CDM 
models the transfer function is: 

-i 

T(k) = {l+[ak + (bk) 1 - 5 + (ck) 2 Y} " (37) 
(Bond e Efstathiou 1984), where a = QA(Q.h 2 )- 1 Mpc) b = 3.0(fi/i 2 )- x Mpc; c = 



1 .7 '(fi/i 2 ) -1 M pc; v — 1.13. It is interesting to note that Eq. (29) is valid only if a « 1, 
since \8\ < 1. This implies than non-linear perturbations, a » 1, must be non-Gaussian. 
In fact when the amplitudes of fluctuations grow up, at a certain point modes are no 
longer independent and start to couple giving rise to non-linear effects that change the 
spectrum and correlation function (Juskiweicz et al 1984 ). There are also some the- 
ories (e.g., cosmic strings (Kibble & Turok 1986)) in which even in the linear regime 
perturbations are not Gaussian. 

1.5 Curvature and isocurvature perturbations. 
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The study of the evolution of density perturbations can be divided into two phases: 

1) perturbations are outside the horizon, in other terms they have a scale A larger than 
Hubble radius r# = ct or A > H^ 1 : 

2) perturbations are inside the horizon, A < H^ 1 . 

In studying the first case it is necessary to use general relativity and one can demonstrate 
that two different kinds of fluctuations exist: curvature or adiabatic (motivated by the 
simplest models of inflation) and isocurvature or isothermal. Curvature fluctuations are 
characterized by a fluctuation in energy density or in space curvature. For them, we can 
write: 

SS 3Sp r Sp m ST dp rn 

-Q = T C - =U l^ 8 ) 

b 4 p r d m 1 p m 

(Efstathiou 1990), where with m has been indicated the matter component, with r 
radiation, with S entropy and with T temperature. Last equation explains why these 
fluctuations are named adiabatic, since for them the entropy variation is zero. 
Isocurvature or isothermal perturbations are not characterized by fluctuations in the 
curvature of the metric, but they are fluctuations in the local equation of state of universe 
and in agreement with the name it results ST = 0. Until fluctuations of isocurvature 
are not inside the horizon, the causality principle does not permit a redistribution 
of energy density. This is possible only when perturbations enter the Horizon and 
isocurvature perturbations can be converted into perturbations in the energy density. 
As a consequence, the distinction between those two kinds of perturbations is no 
longer meaningful after them enter the horizon (Suto et al 1985; Gouda & Sasaki 
1986). The origin of curvature perturbations may be explained inside inflationary 
models or assuming that they are initially present as perturbations of the metric. 
Isocurvature fluctuations may be always produced in inflationary scenarios from fluctu- 
ations in the density number of barions or axions (Efstathiou 1990, Kolb e Turner 1990 ). 

1.6 Perturbations evolution. 

Density perturbations in the components of the universe evolve with time. In order to 
get the evolution equations for S in Newtonian regime, it is possible to use several models. 
In our model, we assume that gravitation dominates on the other interactions and that 
particles (representing galaxies, etc.) move collisionless in the potential <f> of a smooth 
density function (Peebles 1980). The distribution function of particles for position and 
momentum is given by: 

dN = f{x,p,t)d 3 xd 3 p (39) 
and density: 

p(x, t) = ma- 3 J d 3 pf( X , p, t) =p b [l + <S(x, t)] (40) 
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where m is the mass of a particle and pb the background density. Applying Liouville 
theorem to the probability density on a limited region of phase-space of the system we 
have that f verifies the equation: 

| + ^V/-">V^=0 (41) 
at ma* op 

The distribution function f that appears in the previous equations cannot be obtained 
from observations. It is possible to measure moments of f (density, average velocity, etc.). 
We want now to obtain the evolution equations for i5. For this reason, we start integrating 
Eq. ([4l| ) on p and after using Eq. ((icf), we get: 

a3pb di + V J Pfd3p = ° (42) 
If we define velocity as: 

J ma J £_ 

ffdPp 

and introduce it in Eq. (J42) we get: 



(43) 



dS 1 

Pb^r + - V (pv) = (44) 
at a 

We can now multiply Eq. ( f4l| ) for p and integrate it on the momentum: 

4 / Pafd 3 p + — ^ dp [ p a ppfd 3 p + a 3 p(x, t)<t> >a = (45) 
at J ma z J 

this last in Eq. ( ^ ) leaves us with: 

W + 2 ~aft = h V [(1 + 5) V ^ + P~^ d ^ J ™>** P (46) 
and finally using 

< v a v/3 >= r (47) 

mar J Ja d p 

the equation for the evolution of overdensity becomes: 

+ 2 a^ = «? v [(1 + 5) v ^ ] + -* dadp [ (1 + 5) < v ° v " (48) 

(Peebles 1980). The term < v a vp > is the tensor of anisotropy of peculiar velocity. This is 
present in the gradient and then it behaves as a pressure force. If we consider an isolated 
and spherical perturbation, it is possible to assume that initial asymmetries does not 
grow up and so we can suppose, in this hypothesis that < v a vp > = 0. In this case and 
with the linearity assumption 8 « 1 we have: 

^ 

This equation in an Einstein-de Sitter universe (fi = 1, A = 0) has the solutions: 

5+ = A+(x.)ti 5_{x,t) = A-^t- 1 (50) 
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The perturbation is then done of two parts: a growing one, becoming more and more 
important with time, and a decaying one becoming negligible with increasing time, with 
respect to the growing one. 

In the case of open models with no cosmological constant: < 1, A = 0, we can 
write: 

^ = ^Gp(l+(O 1 -l)a), (51) 
and the a(t) evolution can be expressed through the following parametric representation: 



a{r]) = 2(i-°n ) (cosh,? ~ 1) (52) 

t{v) = 2g (i% )3/ 2 (sinh ^ ??) - 

In the case of flat models with positive cosmological constant: < 1, A ^ 0, fl+A/3H^ = 
1, we can write: 

* = ± vG p{i + (nt 1 -i) e P), (53) 

a(t)=(n - 1 -l)- 1/3 sinh 2 /3^yA^ (54) 

Before concluding this section, we want to find an expression for the velocity field in 
the linear regime. Using the equation of motion p = ma 2 x, ^ = — m V <j> and the proper 
velocity of a particle, v = ax, verify the equation: 

— + v- = -— = Gp b a / d d x5(x, t)- (55) 

at a a J |x — x'| 

Supposing that v is a similar solution for the density, v = V+(x, t)t p , we get: 

^ = ^ (56) 
4?r dx a J |x'-x| x ' 

(Peebles 1980). This solution is valid just as that for 5 in the linear regime. At time 

t = to this regime is valid on scales larger than 8h~ 1 Mpc. 



1.7 Non-linear phase. 



Linear evolution is valid only if S « 1 or similarly, if the mass variance, a, is much less 
than unity. When this condition is no longer verified (e.g., if we consider scales smaller 
than 8/i _1 Mpc), it is necessary to develope a non- linear theory. In regions smaller than 
8/i _1 Mpc galaxies are not a Poisson distribution but they tend to cluster. If one wants 
to study the properties of galactic structures or clusters of galaxies, it is necessary to 
introduce a non-linear theory of clustering. A theory of this last item is too complicated to 
be developed in a purely theoretical fashion. The problem can be faced assuming certain 
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approximations that simplifies it (Zel'dovich 1970) or as often it is done, by using N-Body 
simulations of the interesting system. The approximations are often used to furnish the 
initial data to simulations. In the simulations, a large number of particles are randomly 
distributed in a sphere, in the points of a cubic grid, in order to eliminate small scale 
noise. The initial spectrum is obtained perturbing the initial positions by means of a 
superposition of plane waves having random distributed phases and wave vector (West 
et al. 1987). Obviously, the universe is considered in expansion (or comoving coordinates 
are used), and then the equation of motion of particles are numerically solved. For what 
concerns the analytical approximations one of the most used is that of Ze'ldovich (1970). 
This gives a solution to the problem of the grow of perturbations in an universe with 
p = not only in the linear regime but even in the mildly non-linear regime. In this 
approximation, one supposes to have particles with initial position given in Lagrangian 
coordinates q. The positions of particles, at a given time t, are given by: 

x = q + &(<)p(q) (57) 

where x indicates the Eulerian coordinates, p(q) describes the initial density fluctuations 

and b(t) describes their grow in the linear phase and it satisfies the equation: 

d 2 b „ idb da 1 _ , ,„_. 

dfi +2a dtTt=^ Gpb (58) 

The equation of motion of particles, according to the quoted approximation, is given by: 

v = aq + 6p(q) (59) 

The peculiar velocity of particles is given by: 

dx. db , , , i 

while the density of the perturbed system is given by: 

dpk _1 



p(q^) = p 





dqj 




p 


dx k 


= P 



S jk + b(t) _ 

dq 3 



(61) 



Developing the Jacobian present in Eq. ( |6l| ) at first order in 6(t)p(q), one obtains: 

^«-K*)Vqp(q) (62) 

This equation can be re-written, separating the space and time dependence, as in the 
equation for u, and writing: 

b(t) = ti p(q) = V i—^A k exp(ikq) (63) 

k l k l 



in the form: 

Sp = 

k 

(Efstathiou 1990), that leads us back to the linear theory. In other words, Ze'ldovich 
approximation is able to reproduce the linear theory, and is also able to give a good 



= V A k thxp(ikq) (64) 
P ^ 
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approximation in regions with & >> 1. Using the expression for p(q), the Jacobian in 
Eq. ( |6l| ) is a real matrix and symmetric that can be diagonalized. With this p(q) the 
perturbed density can be written as: 

P(q ' t] = (l - &(<)Ai(g))(l - b(t)X 2 ( q ))(l - b(t)X 3 ( q )) (65) 
where Ai, A2, A3 are the three eigenvalues of the Jacobian, describing the expansion and 
contraction of mass along the principal axes. From the structure of the last equation, 
we notice that in regions of high density Eq. ( p5| ) becomes infinite and the structure 
of collapse in a pancake, in a filamentary structure or in a node, according to values 
of eigenvalues. Some N-body simulations (Efstathiou e Silk, 1983) tried to verify 
the prediction of Ze'ldovich approximation, using initial conditions generated using 
a spectrum with a cut-off at low frequencies. The results showed a good agreement 
between theory and simulations, for the initial phases of the evolution (a(t) — 3.6). 
Going on, the approximation is no more valid starting from the time of shell-crossing. 
After shell-crossing, particles does not oscillate any longer around the structure but 
they pass through it making it vanish. This problem has been partly solved supposing 
that particles, before reaching the singularity they sticks the one on the other, due to a 
dissipative term that simulates gravity and then collects on the forming structure. This 
model is known as "adesion-model" (Gurbatov et al. 1985). 

Summarizing, Zel'dovich approximation gives a description of the transition between 
linear and non-linear phase. It is expecially used to get the initial conditions for N-body 
simulations. 



1.8 Quasi-linear regime. 

We have seen in the previous section that in the case of regions of dimension smaller 
than 8h~ 1 Mpc, the linear theory is no more a good approximation and a new theory 
is needed or N-body simulations. Non-linear theory is able to calculate quantities as 
the formation redshift of a given class of objects as galaxies and clusters, the number 
of bound objects having masses larger than a given one, the average virial velocity 
and the correlation function. It is possible to get an estimate of the given quantities as 
that of other not cited, using an intermediate theory between the linear and non-linear 
theory: the quasi-linear theory. This last is obtained adding to the linear theory a 
model of gravitational collapse, just as the spherical collapse model. Important results 
that the theory gives is the bottom- up formation of structures (in the CDM model). 
Other important results are obtained if we identify density peaks in linear regime with 
sites of structure formation. Two important papers in the development of this theory 
are Press-Schechter (1974) and that of Bardeen et al. (1986). This last paper is an 
application of the ideas of the quasi-linear theory to the CDM model. The principles of 
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this approach are the following: 

— Regions of mass larger than M that collapsed can be identified with regions where 
the density contrast evolved according to linear regime, S(M,x), has a value larger 
than a threshold, S c . 

— After collapse regions does not fragment. 

The major drawbacks of the theory, as described in Bardeen et al. (1986) are fundamen- 
tally the fact that the estimates that can be obtained by means of this theory depends 
on the threshold S c , on the ratio between the filtering mass and that of objects and from 
other parameters. Nevertheless, this theory has helped cosmologists in obtaining estimate 
of important quantities as those previously quoted, and at same time give evidences that 
leads to exclude very low values for spectrum normalization. 
1.9 Spherical Collapse 

Spherical symmetry is one of the few cases in which gravitational collapse can be 
solved exactly (Gunn & Gott 1972; Peebles 1980). In fact, as a consequence of Birkhoff's 
theorem, a spherical perturbation evolves as a FRW Universe with density equal to the 
mean density inside the perturbation. 

The simplest spherical perturbation is the top-hat one, i.e. a constant overdensity S 
inside a sphere of radius R; to avoid a feedback reaction on the background model, the 
overdensity has to be surrounded by a spherical underdense shell, such to make the total 
perturbation vanish. The evolution of the radius of the perturbation is then given by a 
Friedmann equation. 

The evolution of a spherical perturbation depends only on its initial overdensity. In an 
Einstein-de Sitter background, any spherical overdensity reaches a singularity (collapse) 
at a final time: 



In an open Universe not any overdensity is going to collapse: the initial density contrast 
has to be such that the total density inside the perturbation overcomes the critical 
density. This can be quantified (not exactly but very accurately) as follows: the growing 
mode saturates at b(t) — 5/2(Qq 1 — 1), so that a perturbation ought to satisfy Si > 
1.69 • 2(f2 ~ 1 — l)/5 to be able collapse. 

Of course, collapse to a singularity is not what really happens in reality. It is typically 
supposed that the structure reaches virial equilibrium at that time. In this case, argu- 
ments based on the virial theorem and on energy conservation show that the structure 




(66) 



By that time its linear density contrast reaches the value: 




(67) 
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reaches a radius equal to half its maximum expansion radius, and a density contrast 
of about 178. In the subsequent evolution the radius and the physical density of the 
virialized structure remains constant, and its density contrast grows with time, as the 
background density decays. Similarly, structures which collapse before are denser than 
the ones which collapse later. 

Spherical collapse is not a realistic description of the formation of real structures; 
however, it has been shown (see Bernardeau 1994 for a rigorous proof or Valageas 2002a, b) 
that high peaks (> 2a) follow spherical collapse, at least in the first phases of their 
evolution. However, a small systematic departure from spherical collapse can change the 
statistical properties of collapse times. 

Spherical collapse can describe the evolution of underdensities. A spherical undcr- 
density is not able to collapse (unless the Universe is closed!), but behaves as an open 
Universe, always expanding unless its borders collide with neighboring regions. At vari- 
ance with overdensities, underdensities tend to be more spherical as they evolve, so that 
the spherical model provides a very good approximation for their evolution. 

1.10 Mass function. 

One of the most important quantities in cosmology is the mass function or multiplicity 
function. It can be described by the relation: 

dN = n(M)dM (68) 

that is the number of objects per unit volume, having a mass in the range M ed M + dM. 
The multiplicity function can also be used to define the luminosity function after having 
fixed the ratio • Obtaining the mass function starting from that of luminosity is com- 
plicated since the ratio ^ is known with noteworthy uncertainty and it is different for 
different objects and moreover the luminosity function for objects like galaxies depend 
on the morphological type (Binggeli & Tamman 1985). Finally trying to determine the 
luminosity function observatively is problematic (see for example G.Efstathiou, R.S.Ellis, 
B.A.Peterson (1988)). 

For the above reasons, the theoretical determination of the mass function is very impor- 
tant. One of the most successful study in the subject is that of Press & Schechter (1974). 
This theory is based upon these hypotheses: 

— The linear density field is described by a stochastic Gaussian field. The statistics of 
the matter distribution is Gaussian. 

— The evolution of density perturbations is that described by the linear theory 
Structures form in those regions where the overdensity linearly evolved and filtered 
with a top-hat filter exceeds a threshold <5 C (6 C — 1.68, obtained from the spherical 
collapse model (Gunn & Gott 1972)). 
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— for 5 > 8 C regions collapse to points. The probability that an object forms at a certain 
point is proportional to the probability that the point is in a region with S > S c given 
by: 



(69) 



f°° 1 f S 2 

P(S, S c ) = / dS L exp - — 
Js c a(2n)2 \ 2o l 

The multiplicity function is given by: 

p{M, z) = - Po ^dM = n(M)MdM (70) 

If we add the conditions £1 = 1, |<5 fe | 2 oc k n , the Press-Schechter solution is autosimilar 
and has the form: 

p ( n + 3\ ( M 

d 4 



1 M n+3 



xexp 

where M*(z) oc (1 + z)~~^ . Several are the problems of the theory: 



— Statistical problems: in the limit of vanishing smoothing radii, or of infinite vari- 
ance, the fraction of collapsed mass, asymptotes to 1/2. This is a signature of linear 
theory: only initially overdense regions, which constitute half of the mass, are able to 
collapse. Nonetheless, underdense regions can be included in larger overdense ones, 
or, more generally, non-collapsed regions have a finite probability of being included 
in larger collapsed ones; this is commonly called cloud-in- cloud problem. PS argued 
that the missing mass would accrete on the formed structures, doubling their mass 
without changing the shape of the MF; however, they did not give a true demonstra- 
tion of that. Then, they multiplied their MF by a "fudge factor" 2. Other authors 
(see Lucchin 1989) used to multiply the MF by a factor (1 + /), with / denoting the 
fraction of mass accreted by the already formed structures. 

— Dynamical problems: the heuristic derivation of the PS MF bypasses all the compli- 
cations related to the highly non-linear dynamics of gravitational collapse. Spherical 
collapse helps in determining the 5 C parameter and in identifying collapsed structures 
with virialized halos. However, the PS procedure completely ignores important dy- 
namical elements, such as the role of tides and the transient filamentary geometry 
of collapsed structures. Moreover, supposing that every structure virializes just after 
collapse is a crude simplification: when a region collapses, all its substructure is sup- 
posed by PS to be erased at once, while in realistic cases the erasure of substructures 
is connected to the two-body interaction of already collapsed clumps, an important 
piece of gravitational dynamics which is completely missed by the PS procedure. 

— Geometrical problems: to estimate the mass function from the fraction of collapsed 
mass at a given scale it is necessary to relate the mass of the formed structure to the 
resolution 
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In practice, the true geometry of the collapsed regions in the Lagrangian space (i.e. as 
mapped in the initial configuration) can be quite complex, especially at intermediate 
and small masses; in this case a different and more sophisticate mass assignment ought 
to be developed, so that geometry is taken into account. For instance, if structures 
are supposed to form in the peaks of the initial field, a different and more geometrical 
way to count collapsed structures could be based on peak abundances. 

Despite all of its problems, the PS procedure proved successful, as compared to N- 
body simulations, and a good starting point for all the subsequent works on the subject 
(Efstathiou et al. (1988); Efstathiou & Rees (1988), Bond ct al. (1991), White, Efstathiou 
& Frenk (1993), Jain & Bertschinger (1994), Lacey & Cole (1994), Efstathiou (1995), 
Bond & Myers (1996b). Most authors reported the PS formula to fit well their N-body 
results; nonetheless, all the authors agree in stating the validity of the PS formula to be 
only statistical, i.e. the existence of the single halos is not well predicted by the linear 
overdensity criterion of PS (see in particular Bond et al. 1991)). There are however some 
exceptions to this general agreement: Brainerd & Villumsen (1992) reported their MF, 
based on a CDM spectrum, to be very similar to a power-law with slope —2, different 
from the PS formula both at small and at large masses. Jain & Bertschinger (1994), 
Gelb & Bertschinger (1994) and Ma & Bertschinger (1994) noted that, to make the PS 
formula agree with their simulations (based on CDM or CHDM spectra), it is necessary 
to lower the value of the S c parameter as redshift increases. The same thing was found by 
Klypin et al. (1995), but was interpreted as an artifact of their clump-finding algorithm. 
Recent simulations seem to confirm this trend. 

Lacey & Cole (1994) extended the comparison to N-body simulations to the pre- 
dictions for merging histories of dark-matter halos; they found again a good agreement 
between theory and simulations. This fact is noteworthy, as merging histories contain 
much more detailed information about hierarchical collapse. Several improvements of 
the theory exists: Bond et al. 1991; Lacey & Cole (1993) (the so called Extended Press 
& Schechter (EPS) formalism; Sheth & Tormcn (2002). 

1.11 CDM and HDM cosmogony. 

The study of origin and formation of structure in the universe has been historically 
fundamentally framed into two theories: the CDM theory, in which WIMPS constitutes 
the main part of Dark Matter, and HDM in which neutrinos dominate. As we are going 
to see, structure formation in these scenarios is completely different since WIMPS and 
neutrinos are subject to different physical phenomena and then the transfer function is 
noteworthy different in the two cases. Both theories have the same starting points: 
1) The universe is fundamentally constituted by Dark Matter (WIMPS in the CDM 
model and neutrinos in the HDM) and p = p c ( O = 1). 
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2) Baryons give a small contribution to the mass of the universe. 

3) Fluctuations originating in the primordial universe are adiabatic, scale invariant, 
n = l, and Gaussian (Efstathiou 1990). 

If universe is dominated by neutrinos with mass m„ = 30ew the transfer function is 
determined by the free-streaming (or Landau damping) of neutrinos. This phenomenon 
consists in the smoothing of inhomogencities in the primordial universe (due to perturba- 
tions in the acollisional components) because of the motion of neutrinos from overdense 
to underdense regions. Neutrinos diffusion and the smoothing of inhomogeneities is possi- 
ble only before t = t eq . After this epoch, there is no longer free-streaming but the density 
perturbation has definitely changed by its previous action. Free-streaming scale or mass 
can be estimated calculating the distance covered from a particle decoupled from plasma. 
Results that one obtains for the free-streaming scale and mass is (Kolb & Turner 1990): 

X FS w AOMpcim^/SOev)- 1 M FS w 10 15 (m ly /30ew)" 2 M (72) 

Because of free-streaming of neutrinos the HDM spectrum is characterized by a cut-off 
at short wave-length. The result is that the first objects that form are superclusters and 
structure formation proceeds because of fragmentation. Zcl'dovich (1970) showed that 
the first structure to form are flat and were called " pancakes" . After formation, these 
objects enter in the non-linear phase along one of the axes and baryons inside start to col- 
lide and dissipate their gravitational energy. Galaxies form for fragmentation processes. 
Structure formation follows a 'top-down' scheme, that means that larger objects (e.g., 
clusters) form before, while smaller objects (e.g., galaxies) later. N-body simulations of 
HDM universes (Klypin c Shandarin 1983; Frenk et al 1983; Ccntrclla and Mellot 1983) 
showed that on scales larger than 10 Mpc structures is qualitatively similar to voids and 
to the filamentary structure that is visible in the CFA, but the clustering measured in 
N-body simulations is larger than that observed in the CFA. When one tries to reproduce 
the observed correlation function, one arrives to the conclusion that pancaking should 
have happened at redshift z < 1, in disagreement with observed galaxies having z > 1 
and QSO with z > 3. A further problem of the model is that of the peculiar velocity that 
are smaller than values obtained from observations (Kolb e Turner 1990). 
The HDM model after a series of studies in the '80s has been abandoned for the problems 
it has and replaced by another model, the CDM, which is in better agreement with ob- 
servations. The CDM model has a spectrum without a cut-off at short wave -length (at 
least till scales much smaller than galactic scales) because the damping scale is unimpor- 
tant for WIMPS with mass > 1 Gev. Structure formation is typically hyerachical: from 
smaller scale structure to larger ones. This scheme is a 'bottom-up' scheme. When the 
CDM model was introduced it obtained noteworthy successes in the description of the 
characteristic of the universe (clustering statistics of galaxies, peculiar velocities, CMBR 
fluctuations) from the galactic scale on ( Peebles 1982 b; Blumenthal et al 1984; Bardeen 
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et al 1986; White et al 1987; Frenk et al 1988; Efstathiou 1990). The model has shown 
some weak points, when compared with more and more precise data. 
The reason of the success of the CDM model is fundamentally due to the fact that 
WIMPS interact with matter by means of gravity only, and does not feel the effect of 
pressure forces due to interaction with radiation (to which matter components are sub- 
ject). Structure formation starts before in the CDM component, at t < t eq , which give 
rise to the potential wells in which baryonic matter can then fall. It is important to notice 
that, in order to reproduce observations, an additive hypothesis must be added: the bias- 
ing hypothesis, that can be summarized in: light does not trace mass. Q Typical problems 
of the CDM model in absence of biasing are the too high values for the r.m.s. of peculiar 
velocity of couple of galaxies (values of 1000 km/s vs. 300± 50 km/s observed). Another 
problem is the correlation length, r$, in N-body simulations, for the correlation function 
£ which is equal to 1.3ft, _1 Mpc, smaller than the observed value 5.5/i _1 Mpc ( Davis et al. 
1985). On the other side, if, in order to eliminate the quoted problems one introduce the 
biasing hypothesis, there is the supplementary problem of finding a physical mechanism 
that explains the origin of bias. Several conjectures has been proposed (Rees 1985; Silk 
1985; Dekel e Rees 1987) but there is no full agreement on them. 

Summarizing, one can tell that although at the beginning the standard form of 
CDM was very successful in describing the structures observed in the Universe (galaxy 
clustering statistics, structure formation epochs, peculiar velocity flows) (Peebles, 1982; 
Blumenthal et al. 1984; Bardeen et al. 1986 -hereafter BBKS; White et al. 1987; Frenk 
et al. 1988; Efstathiou 1990) recent measurements have shown several deficiencies in the 
model, at least, when any bias of the distribution of galaxies relative to the mass is con- 
stant with scale (see Babul & White 1991; Bower et al. 1993; Del Popolo & Gambera 
1998a, b). Some of the most difficult problems that must be reconciled with the theory 
are: 

— the magnitude of the dipole of the angular distribution of optically selected galaxies 
(Lahav et al. 1988; Kaiser & Lahav 1989); 

— the possible observations of clusters of galaxies with high velocity dispersion at z > 0.5 
(Evrard 1989); 

— the strong clustering of rich clusters of galaxies, £, C c(r) — (r/25h~ 1 Mpc)~ 2 , far in 
excess of CDM predictions (Bahcall & Soneira 1983); 

— the X-ray temperature distribution function of clusters, over-producing the observed 
clusters abundances (Bartlett & Silk 1993); 

— the conflict between the normalization of the spectrum of the perturbation which is 
required by different types of observations; 

5 If this hypothesis is not introduced the discrepancies between data and theory can be reduced 
assuming a value of Hq « 25km /M pes 
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— the incorrect scale dependence of the galaxy correlation function, £(r), on scales 10 to 
100 /i _1 Mpc, having £(r) too little power on the large scales compared to the power on 
smaller scales (Lahav et al. 1989; Maddox et al. 1990a; Saunders et al. 1991; Peacock 
1991; Peacock & Nicholson 1991). 

— Normalization obtained from COBE data (Smoot et al. 1992) on scales of the order 
of 10 3 Mpc requires as = 0.95 ± 0.2, where as is the rms value of in a sphere of 
8/i _1 Mpc. Normalization on scales 10 to 50Mpc obtained from QDOT (Kaiser et al. 
1991) and POTENT (Dekel et al. 1992) requires that a 8 is in the range 0.7 -j- 1.1, 
which is compatible with COBE normalization while the observations of the pairwise 
velocity dispersion of galaxies on scales r < 3Mpc seem to require as < 0.5. 

— Another problem of CDM model is the incorrect scale dependence of the galaxy 
correlation function, £(r), on scales 10 to 100 Mpc, having £(r) too little power on 
the large scales compared to the power on smaller scales. The APM survey (Maddox 
et al. 1990), giving the galaxy angular correlation function, the 1.2 Jy IRAS power 
spectrum, the QDOT survey (Saunders et al. 1991), X-ray observations (Lahav et al. 
1989) and radio observations (Peacock 1991; Peacock & Nicholson 1991) agree with 
the quoted conclusion. As shown in recent studies of galaxy clustering on large scales 
(Maddox et al. 1990; Efstathiou et al. 1990b; Saunders et al. 1991) the measured 
rms fluctuations within spheres of radius 20h~ 1 Mpc have value 2-3 times larger than 
that predicted by the CDM model. 

— Density profiles of CDM halos: the cusp obtained from numerical simulations seems 
too steep. - 

— simulations might yield too many satellites for galaxies like our own. Though this 
second problem may have been the result of bad comparison of simulations with 
observations. This yielded a surge of interest in the last 2-3 years for Warm Dark 
matter and slightly collisional matter. 

These discrepancies between the theoretical predictions of the CDM model and 
the observations led many authors to conclude that the shape of the CDM spectrum 
is incorrect and to search alternative models (Peebles 1984; Shafi & Stecker 1984; 
Valdarnini & Bonomctto 1985; Holtzman 1989; Efstathiou et al. 1990a; Turner 1991; 
Schaefer 1991; Cen et al. 1992; Schaefer & Shafi 1993; Holtzman & Primack 1993; Bower 
et al. 1993). 

Alternative models with more large-scale power than CDM have been introduced in 
order to solve the latter problem. Several authors (Peebles 1984; Efstathiou et al. 1990a; 
Turner 1991) have lowered the matter density under the critical value (Cl m < 1) and 
have added a cosmological constant in order to retain a flat Universe (Q m + Cl\ = 1) . 
The spectrum of the matter density is specified by the transfer function, but its shape is 
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affected because of the fact that the epoch of matter-radiation equality is earlier, 1 + z eq 
being increased by a factor l/f2 m . The epoch of matter-radiation equality is earlier, 
because 1 + z eq is increased by a factor 1/Q m . Around the epoch z\ the growth of the 
density contrast slows down and ceases after z\. As a consequence the normalization 
of the transfer function begins to fall, even if its shape is retained (and pushes its 
imprint to larger scales). Mixed Dark Matter models (MDM) ( Bond et al. 1988; Shah 
& Stecker 1984; Valdarnini & Bonometto 1985; Schaefer et al. 1989; Holtzman 1989; 
Schaefer 1991; Shaefer & Shafi 1993; Holtzman & Primack 1993) increase the large-scale 
power because neutrinos free-streaming damps the power on small scales. Alternatively, 
changing the primeval spectrum several problems of CDM are solved (Cen et al. 1992). 
Finally it is possible to assume that the threshold for galaxy formation is not spatially 
invariant but weakly modulated (2% — 3% on scales r > \0h~ 1 Mpc) by large scale 
density fluctuations, with the result that the clustering on large-scale is significantly 
increased (Bower et al. 1993). 

Conclusions. 

This paper provides a review of the variants of dark matter (CDM, HDM) which 
are thought to be fundamental components of the universe and their role in origin and 
evolution of structures. Studies of several decades have shown that, if we have a right 
knowledge of the law of gravity, dark matter is a fundamental component of our universe. 
While models based upon Hot Dark Matter (e.g., neutrinos) gives a reasonable description 
of structures on large scales models based upon Cold Dark Matter (e.g., axions) are 
more successful in describing small and intermediate scales. A fundamental ingredient 
in the recipe of structure formation is inflation which provides a spectrum of adiabatic 
Gaussian perturbations which can be well described by a power-law spectrum, tilted 
from the Harrison-Zel'dovich spectrum, normally tilted so as to provide extra large scale 
power. The magnitude of the tilt may be modest or pronounced. The details of structure 
formation are very sensitive to the matter content of the universe. It appears that if cold 
dark matter is the main constituent of the universe, present observations require that 
the initial perturbations be adiabatic isocurvaturc perturbations generate excessively 
large cmb anisotropies for the same final density perturbation. Adiabatic perturbations 
are exactly what inflation provides. In CDM models, the only remaining alternative would 
appear to be texture seeded models, which have been placed in jeopardy by a combination 
of microwave anisotropy and velocity data, though the death blow apparently remains to 
be struck (Pen, Spergel & Turok 1992). The survey was completed by examining variants 
on the CDM model which may be better suited to explaining the observational data. The 
standard technique is to utilize additional matter (be it a component of hot dark matter 
or of a cosmological constant) to remove short-scale power from the CDM spectrum. Hot 
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dark matter does this by free-streaming, a cosmological constant by delaying matter- 
radiation equality. Because this power can be removed over a much shorter range of 
scales than with tilt, it enables an explanation of the observed deficit of short-scale 
power relative to intermediate scale power in the spectrum. 

MDM (Mixed dark Model) adds yet another new parameter, roughly speaking an abil- 
ity to remove short-scale power from the spectrum while leaving large scales untouched, 
and may be necessary should all present observations stand up. It appears likely that 
MDM will however need an initial spectrum close to n = 1 with no gravitational waves 
if it is to succeed. Some studies (Peebles 1984; Efstathiou et al. 1990a; Turner 1991) has 
shown that lowering the matter density under the critical value (Cl m < 1) and adding 
a cosmological constant in order to retain a flat Universe (f2 m + Qj\ = 1), gives good 
results in the case of Q m = 0.3. Moreover new observational evidences (see Corasaniti, 
Copeland 2001) indicates that we are living in a A ^ universe. 
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